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T. E. Harris was a pioneer par excellence in many fields of proba- 
bility theory. In this paper, we give a brief survey of the many funda- 
mental contributions of Harris to the theory of branching processes, 
starting with his doctoral work at Princeton in the late forties and 
culminating in his fundamental book "The Theory of Branching Pro- 
cesses," published in 1963. 

1. Introduction. T. E. Harris wrote the first definitive book [5] on branch- 
ing processes, published in 1963. It covered much of the work on the subject 
up to that time, a sizeable part due to Harris himself. It identified the subject 
of branching processes and resulted in a great deal of interest in the subject, 
among both mathematicians and statisticians. Between 1963 and 1970, a 
vast number of papers on branching processes appeared in many good jour- 
nals specializing in probability theory and mathematical statistics, and by 
1971 more books on the subject appeared both in the U.S. and elsewhere 
[15, 18]. Harris himself moved on to work on other beautiful topics such as 
percolation and interacting particle systems. As with branching processes, 
his work in these other areas was profound. T. E. Harris was pioneer par 
excellence, creating many areas of research in which he laid the foundations 
that others built on. In what follows, we present a brief account of Harris's 
contribution to branching processes. 

Harris's 1947 PhD dissertation at the mathematics department of Prince- 
ton University was on branching processes, titled "Some theorems on Bernoulli 
multiplicative processes." This was followed in 1948 by his basic paper [6] 
in the Annals of Mathematical Statistics. In [6], he used the term branching 
processes, a term which had also been used by Russian mathematicians; he 
treated the single type discrete time branching process. He also coined the 
term Galton-Watson branching process for this process. His main focus in 
[6] was on the supercritical case; we now give a description of this work. 
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2. Single type, discrete time case. Let {pj}j>o be a probability distri- 
bution. Let {£,n,k]n > 0,/c > 1} be an array of nonnegative integer valued 
random variables that are i.i.d. (independent and identically distributed) 
with distribution {pj}j>o- Let Zq be a positive integer. Now set 

Z 

(1) ^i = J>,fc 

k=l 

and for n > 1, Z n+ \ = J2k=i £n,k if > an d if Z n = 0. Then the sequence 
{Z n } n >o is called a Galton- Watson branching process with initial population 
Zq and offspring distribution {pj}j>o- Clearly, {Z n } n >Q is a Markov chain 
with time homogeneous transition probabilities and the nonnegative integers 
as the state space. The transition probabilities are given by 

Pij = P ^ r = J J for i > 1 and p o = 1, 

where {£ r }r>i are i.i.d. with distribution {pk}k>o- 

One can interpret the sequence {Z n } n >Q as follows. If Z n is thought of as 
the number of individuals in the nth generation, then each one of them pro- 
duces a random number of children with distribution {pj}j>o independently 
of others in the nth generation as well as any past ancestors. The total 
number Z n+ \ of all these individuals is the size of the (n + l)st generation. 

An important parameter in determining how the sequence {Z n } n >Q be- 
haves for n large is the offspring mean m = Y2j JPj ■ Here are some basic 
results. 

Theorem 2.1. Let < m = Y.JLiJPj < 00 an d ^ et P(® < z o < °°) = !• 
Then: 

(i) m < 1 P(Z n — > as n — > oo) = 1, 

(ii) m = l,pi < 1 => P(Z n — >• as n — > oo) = 1, 

(hi) m > 1 => P(^ n — >• as n — >■ oo | Zq = 1) = g < 1, 

where q is the unique root of the equation 

oo 

s = /0) = X>i sj > 0<s<l. 

j=0 

Further, P(Z n — > oo as n — > oo | Zq = 1) = 1 — q, and for any k > 1, P(Z n — » 
as n — )• oo | = A;) = g fc . 

Harris in his book [5] notes in that in 1874 in [19], Galton and Watson 
did notice that the extinction probability q satisfied q = f(q), but failed to 
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notice that if m > 1 the relevant root is less than one. Galton and Watson's 
work was motivated by the problem of the survival of British peerage names, 
posed by Galton in the London Times in the 1870s. 

In his paper [6], which is based on his doctoral thesis, Harris focused 
mainly on the supercritical case, that is, m > 1. The case m = 1 is called 
critical and m < 1 is the subcritical case. Let {pj},m,{Z n } be as in Theo- 
rem 1. 

Theorem 2.2 (Supercritical case, [6, 7]). Assume po = 0,pi < l,m > 1, 
Y^=ij 2 Pj < 00 an d < < oo- £ef W n = Z n /m n ,n > 0. TTien there exists 
a nonnegative random variable W such that: 

(i) E((W n - W) 2 | Z ) -4 as n -> oo, 

(ii) P(V = 0)=0, 

(hi) TV /ias an absolutely continuous distribution an (0, oo) mi/t a continu- 
ous density, 
(iv) S(W|Z = 1) = 1. 

Harris [5] observes that J. L. Doob seems to have been the hrst to note that 
{H^„} n >o is a martingale and, being nonnegative, converges a.s. as n — > oo. 
Kesten and Stigum [13] improved on this, as follows. 

Theorem 2.3 [13]. Let p = 0,pi < 1,0 < Z < oo, 1 < m, and W n = 
Then: 

(i) J*C°Sj)Pj < °o =^ W n — t- W a.s. and in mean, where P(W = 0) = 
0, E(W | Zq = 1) = 1 and W has an absolutely continuous distribution on 
(0,oo). 

(ii) Ei°i( lo gi)Pi = °°=> W n ,->0, a.s. 

The work of A. N. Kolmogorov [14] in 1938 and A. M. Yaglom [20] in 
1947 (see [5]) led to the following. 

Theorem 2.4 (Critical case). Suppose m = \,p\ < 1 and Y1T j 2 Pj < 00 • 
Then, as n^oo, 

(i) nP(Z n >0\Z = 1)^ 2 }, where a 2 = E? fPj ~ h 

(ii) > x | Z = 1, Z n > 0) -> e~ 2 / CT2:l ', /or a22 < x < oo. 

Theorem 2.5 (Subcritical case). Lei ?n < 1. Then for all j > 1, lim n P{Z n — 
j\Z n >0) = bj exists, < bj < oo and X^li = !• 

In his book [5], Harris presents extensions of Theorems 2.1, 2.2, 2.4 and 
2.5 to the multitype (finite type) case. In [13], Kesten and Stigum established 
the analog of Theorem 2.3 above for the multitype Galton- Watson process. 
See Athreya and Ney [1] for details; see also Sevastyanov [18] and Mode [15]. 
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3. Single type, age dependent case. In 1948 Harris, with Richard Bell- 
man [3, 10], formulated the theory of age dependent branching processes, 
where each individual lives a random length of time and on death creates a 
random number of individuals, and all individuals live and reproduce inde- 
pendently of each other. Assuming all moments on the offspring distribution 
and an absolutely continuous life time distribution, they established an inte- 
gral equation for the probability generating function of Z(t), the population 
size at time t. They showed that in the supercritical case, Z(t)e~ at converges 
in probability to a limit random variable W, where a is the Malthusian pa- 
rameter defined by m J e~ au dG(u) = 1, with (?(-) being the distribution 
function of the lifetime of an individual. They further showed that W is 
nontrivial and has an absolutely continuous distribution on (0,oo). There 
are analogs of Theorems 2.3 and 2.4 for this CclSG, clS well. 

Conditions for the supercritical case were relaxed by later authors; see 
Athreya and Ney [1]. 

4. General type case. Harris also considered branching processes with 
arbitrary type space by using the point process approach. Here in any gen- 
eration, one has a finite point process an some type space X. The basic 
branching property of independent production is retained. An individual lo- 
cated at x G X produces children according to a point process over X whose 
distribution depends on x. All individuals act independently of each other. 
For this, Harris used the method of moment generating functions. In [9], he 
established the analog of Theorem 2.2 in this context, and applied this to 
nuclear cascades and related processes, as well as a one-dimensional neutron 
model. For details on this, see Harris's book [5]. Harris mentions that J. E. 
Moyal worked on similar ideas. In the 1970s, Jagers and his colloborators in 
Sweden developed this topic further in great detail (see [11]). See also Ney 
[16, 17]. 

5. Cosmic-ray cascades. Harris studied the theory of cosmic-rays cas- 
cades and supplemented the work of nuclear physicists; Chapter 7 of his 
1963 book [5] deals with this topic. We present a brief summary of Harris's 
work on cosmic-ray cascades as discussed in his paper [8] . Here are the model 
assumptions: 

(1) A photon of positive energy e, moving through homogeneous material, 
has probability A dt + o(dt) of being transformed in the thickness interval 
(t,t + dt) into two electrons, positive or negative, which receive energies eU 
and s(l — U), respectively, where U is a random variable with an absolutely 
continuous distribution in (0,1). Note that the role of time parameter is 
played by the thickness of the material. 

(2) An electron loses (by "collision" or "ionization") a deterministic amount 
of energy (3t in an interval of length t. 
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(3) An electron radiates photons in such a way that the probability that 
an electron of energy e emits a photon of energy between eu and e(u + du) in 
a small thickness interval of length dt is k(u) dudt. Further, the energy that 
goes to the radiated photon is subtracted from that of the parent photon. 
A special case of interest for k(-) is k(u) = ^ + ko(u) with | I < c (l — 
u)~ h ,0 < u < 1, where c and [i are constants and b < 2. 

Here A and (3 are constants independent of t and e. 

Under the above assumptions, Harris shows that if f3 = and £o(t) is 
the energy at time t of an electron with eo(0) = 1 then, for t > 0,Xo(t) = 
— log £o (t) nas an infinitely divisible distribution with probability density 
ht(x),x > 0, whose characteristic function (Fourier transform) is given by 

POO / POD \ 

J e i9x h t (x)dx = exp(t J (e i0u - l)k(l - e- u )e~ u du) . 

Harris observes that the special case when k(u) = — log ^„ tt ^ an d fh{x) is the 

Gamma density — > this was given by Bhabha and Heitler [4]. 

Next, Harris considers the random process N(e,t), the total number of 
electrons at t whose energies are greater than e, for e > 0. Let fi(s,e,t) = 
Ei(s N ^ £ '^),0 < s < 1, where Ei stands for expectation when the starting 
particle is of energy 1, and is a photon for % = 1 and an electron for i = 2. 
Harris shows that the following integro-differential equation holds: 



|(, £l «)= J [ , [/ I (,V)/ ! (» 1 ^, 1 ) 



f2(s,e,t) 



k{u) du, 



with fi(s,e, 0) = /i(s,l,i) = f2(s,l,t) = 1 for t > and /2(s,e, 0) = s for 
e < 1. Harris shows that the earlier results of Bartlett and Kendall [2] and 
of Janossy [12] could be deduced from the above. 

Harris introduces a vector valued Markov process (I(t),C(t)),t > 0, where 
I(t) is the condition of a single particle at time t which can be a photon 
(/ = 1) or an electron (/ = 2) and has energy C(t). He then derives the 
limiting distribution of the process (I(t)X{t)) as t — > oo (assuming j3 = 0) 
and is able to deduce the earlier results of other authors as special cases. 

Harris also obtained results for cascades with f3 > 0. In particular, he 
shows that when (3 > 0, the energy £\{t) of an electron at time t can be 
represented by 



E 1 (t)=max{o,e (t)(l-p J 



6. Concluding remarks. T. E. Harris was deeply involved in the devel- 
opment of all aspects of contemporary branching process theory. He laid a 
rigorous foundation to areas where it had been lacking. His 1963 book [5] 
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is a beautiful and major work of scholarship. One can substantially credit 
to its publication the explosion of work on branching processes in the 1960s 
and 1970s and up to the present. It set the impetus and direction of research 
on the subject for many years. The present authors owe T. E. Harris a deep 
debt of gratitude for this. 
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